Symmetry relations in chemical kinetics arising from microscopic reversibility 
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It is shown that the kinetics of time-reversible chemical reactions having the same equilibrium 
constant but different initial conditions are closely related to one another by a directly measurable 
symmetry relation analogous to chemical detailed balance. In contrast to detailed balance, however, 
this relation does not require knowledge of the elementary steps that underlie the reaction, and 
remains valid in regimes where the concept of rate constants is ill-defined, such as at very short times 
and in the presence of low activation barriers. Numerical simulations of a model of isomerization in 
solution are provided to illustrate the symmetry under such conditions, and potential applications 
in protein folding-unfolding are pointed out. 



Consider a reversible unimolecular reaction of the type 



X 



A ^ B, 



(1) 



either in solution or in the presence of a nonreactive 
buffer gas, at low concentrations of species A and B. 
When the above reaction is elementary, i.e. when it pro- 
ceeds without the formation of intermediate complexes, 
the so-called detailed balance relation establishes that the 
kinetic rate constants corresponding to the forward and 
reverse directions of Eq. QJ, kf and k r respectively, are 
related to each other by a constant factor, namely (sec 
e.g. 0) 



(2) 



where K = cb /ca is the "equilibrium constant" given by 
the ratio of the equilibrium concentrations of B and A, 
cb and ca respectively. When intermediate complexes 
are involved, e.g. when the elementary steps are given 
by A <=> I <=> B, where I is an intermediate complex pos- 
sibly involving solvent molecules, the detailed balance re- 
lation still applies to each elementary step Although 
this relation is a tool of proven usefulness in the analy- 
sis of various chemical reactions, its applicability is lim- 
ited to cases where the elementary steps of the reaction 
can be identified, and hinges upon the existence of rate 
"constants" ; during transient times or in the presence of 
low activation barriers, for example, reactions are not ex- 
pected to proceed in "local equilibrium" Q , and thus in 
these regimes rate constants are ill-defined concepts. 

In the present contribution, the assumption of micro- 
scopic reversibility will be exploited to derive an equal- 
ity relating time-dependent concentrations to equilib- 
rium constants in reversible reactions, such as Eq. IjTjI. 
where the time-dependent concentrations are measured 
in two complementary experiments involving different 
initial conditions (see below) . This relation remains valid 
in the aforementioned regimes where rate constants are 
not well-defined, and is independent of the elementary 



steps that underlie the reaction. In the specific context of 
isomerization reactions expressed by Eq. JQ , the equality 
is given by 



,(*) 



c r A (t) 



(3) 



and is satisfied at all times t > after the reaction has 
started at t = 0. In the above equation and subsequently 
in this paper, the superscripts / and r specify two dif- 
ferent experiments: one that starts with a pure sample 
of species corresponding to the left side of the chemi- 
cal equation (in this case, A), and another with a pure 
sample of species corresponding to the right side of the 
equality (£?), respectively. (Note that this definition of 
directionality is distinct from that associated with rate 
constants). Therefore, according to Eq. |J3}, if one ob- 
server starts with a sample containing A only and mea- 
sures the fractional concentration of B as a function of 
time (c B (t)), another observer that starts with a sample 
of B only will observe exactly the same time-dependence 
of the fractional concentration of species A (c r A (t)), up 
to an overall constant factor K — cb/c a (see Fig. ^ for 
a numerical example). This equality can thus be seen as 
a quantitative, nonequilibrium formulation of the princi- 
ple of microscopic reversibility 0, 0] that is particularly 
suited for chemical reactions. 

The result expressed by Eq. © can also be regarded as 
part of a recent endeavor to establish directly measurable 
equalities that relate opposite "directions" of nonequi- 
librium processes in time-reversible systems pL These 
fluctuation theorems have been derived (e.g. " @, @) 
and experimentally verified (e.g. Il0|) in a variety of 
physical contexts, and are set apart from other nonequi- 
librium results as they do not rely on approximate ar- 
guments that are traditionally invoked in the study of 
out-of-equilibrium problems, such as near-equilibrium as- 
sumptions, or - in the present context of chemical kinet- 
ics - transition state theory. 

Before providing the derivation of Eq. ©, note that 
this equality is consistent with known limits and phc- 
nomenological rate equations. Indeed, as i oo, the 
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time-dependent concentrations c B (t) and c r A (t) approach 
their asymptotic equilibrium values, c B and ca respec- 
tively, and hence the ratio approaches the correct equi- 
librium constant in this limit. The limit t — > involves 
the ratio of two vanishing quantities, as the initial con- 
centration of products is zero in both reaction directions. 
However, taking the time-derivative of Eq. 10 at t = 
and introducing the rate constants as 

c f B (0) = k f c{(0) 
c r A (0) = k r c r B (0), 

one obtains the statement of detailed balance, Eq. J5J), 
since c A (0) = c r B (0) = 1 before the reaction takes place 
(recall that concentrations are fractional in this paper). If 
the reaction follows a first-order kinetics (see e.g. Q), it 
can be easily verified that the time-dependent concentra- 
tions satisfy Eq. ©. It is important to emphasize, how- 
ever, that this equality remains valid in regimes where 
these phenomenological descriptions break down, such 
as at short (transient) times or low activation barriers 
(i.e. high temperatures); this property will be illustrated 
through explicit numerical experiments below. 

Two proofs of Eq. J2J will now be offered. The first 
shows that this result follows almost directly from the 
property of microscopic reversibility (according to On- 
sager |ll|; see also 0,0]), which states that - in equilib- 
rium - the joint probability of observing a function /(q) 
of some configurational degrees of freedom q of the sys- 
tem taking the value /(q) = A at time t = and the 
value /(q) = B at a later time t — r is equal to the 
probability of observing the reverse, namely /(q) = B at 
time t = and /(q) = A at time t = r. Mathematically, 
this property can be summarized as 

Peq (A,0-B,T) = p eq (B,0;A,r), 

in obvious notation. In particular, /(q) could be a func- 
tion that indicates the chemical species corresponding 
to the state q. From elementary probability theory, we 
can decompose the above joint probabilities in terms of 
marginal and conditional probabilities, so that 

Peq (A)p eq (B, t\A, 0) = Peq (B) Peq (A, T \B, 0), 

a relation that is sometimes known as "detailed balance" 
too (note, however, that this result is more general than 
the "chemical" detailed balance relation, Eq. J2Jl). The 
crucial step to obtain Eq. J3J is to realize that the condi- 
tional probabilities p eq (B, t\A, 0) and p e q(A, t\B, 0) dic- 
tating transitions between A and B in equilibrium coin- 
cide with the nonequilibrium concentrations c b (t) and 
c t a (t) introduced earlier in the paper, respectively, pro- 
vided the initial statistical state of these nonequilibrium 
processes are in local equilibrium (see below); the central 
result of this work then follows by using this identifica- 
tion, and noticing that K = p eq (B)/p eq (A) (see also 0] 



for a similar result and derivation applied in a different 
context). 

The second proof uses the full microscopic description 
of the dynamics, thus revealing the underlying dynami- 
cal and statistical assumptions that lead to Eq. and 
avoiding the last and perhaps unfamiliar step of the first 
proof that identifies equilibrium transition probabilities 
with nonequilibrium concentrations. This derivation as- 
sumes that the dynamics of the system (molecule and en- 
vironment) is Hamiltonian, and for simplicity of notation 
the reaction coordinate is assumed to be one-dimensional. 
The microscopic state of the system is fully specified by 
the vector x = (q, Q,p, P), where q and p are the posi- 
tion and momentum associated with the reaction coor- 
dinate, and Q and P are the remaining configurational 
and kinetic degrees of freedom of the molecule and envi- 
ronment. Chemical species are defined by the values of 
q: Species A corresponds to q < and species B corre- 
sponds to q > 0. 

The initial statistical state of the reaction in either 
direction is assumed to be in "local equilibrium" Q . For 
example, the initial statistical state corresponding to the 
forward process in which no B is present at t = is given 
by 

j _ p eq (H(x))0(-q) 
PoW- Ca 

where ca = fdxp eq (H(x))9(—q) normalizes the distri- 
bution and coincides with the equilibrium concentration 
of A, 9{q) is the step function, and p eg (£f(x)) is the equi- 
librium distribution which is assumed to depend on the 
total Hamiltonian H (x) only (e.g. the canonical or micro- 
canonical distributions) . The case of the reverse reaction 
follows by analogy. Note that this assumption is funda- 
mentally different from the condition of local equilibrium 
at all times t > that underpins transition state theory 
[2|, and in general this initial distribution will evolve to- 
wards states that violate the local equilibrium condition 
at later times. The time-dependent concentration of B 
in the forward process can be written as 

4(t ) =^xo ^ (J?(Xo))g( - 9 ° ) %0, (4) 

J CA 

where subscripts in coordinates indicate time (e.g., qt 
is the q-component of the vector x t obtained by inte- 
grating the equations of motion t units of time with ini- 
tial condition xq). The assumption that the dynamics is 
Hamiltonian excludes the possibility that the coordinates 
x exchange energy with other particles in a larger envi- 
ronment during the reaction time; if the presence of other 
degrees of freedom invalidate this assumption, these can 
be included in Q and P as necessary until one obtains a 
closed system (alternatively, one can adopt a reversible 
Markovian model to mimic the behavior of a system cou- 
pled to a large heat reservoir Jl2| ) . This ensures not only 
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that energy is conserved, i.e. fl(xo) = H(x t ), but also 
that the dynamics preserves the measure in phase space 
[l3|. i.e. dxQ = dx t , so that Eq. can be written as 

4(t)= [d ^ H ^ ))6{ - qo) e{ qt ). (5) 

J C A 

The next step invokes the assumption of dynami- 
cal reversibility which corresponds to the property 
fl(x) = if (x*), where x* = (g, Q, -p, -P) is the "time- 
reversed" counterpart of x, obtained through a simple 
velocity-reversal operation. This property rules out the 
presence of magnetic fields or dissipative terms in the 
Hamiltonian, but is otherwise generally satisfied. It im- 
plies, in particular, that for any trajectory xq • • • x t that 
is a solution of Hamilton's equations of motion, there 
exists another solution Xo • • ■ Xj that coincides with the 
reverse trajectory x£ • • • Xq (see Fig. 1 of Ref. for an 
illustration of this statement). Thus, since the velocity- 
reversal operation preserves the measure in phase space, 
i.e. dx = dx* , one can change the integration variables 
in Eq. (J5J as 

f f p eq {H{x a ))6{-q t ) 

c J B (t) = / dx a % ). 

J CA 

The above integral can be recognized as the definition of 
c r A {t) up to the normalization factor ca/cb] taking this 
into account, one obtains 

4(*) = ^^(t), 

CA 

which is the desired result, Eq. 

In order to illustrate the validity of this result, con- 
sider a one-dimensional isomerization model where the 
reaction coordinate is coupled to a chain of oscillators 
representing a surrounding liquid or buffer gas. The full 
Hamiltonian of the model is given by 

fl(x) = H mo i{p, q) + H env (P, Q) + ff int (g, Q), 

where 

2 i 

H mo l(p, q) = y + a(q 2 - l) 2 + -(q - l) 2 

is the molecule Hamiltonian with a double well poten- 
tial whose barrier height and degree of asymmetry are 
controlled by the parameters a and b respectively (6 = 
corresponds to a symmetric double well), 

fle^P, Q) = ^ + E T + E ^ "f " l)4 
i=l i=2 

is the Hamiltonian of the environment, and 

N 

fl»t(g,Q)=ffg 4 E^ 

i=l 




Time (t) 

FIG. 1: Numerical results for the one-dimensional isomeriza- 
tion model described in the text, illustrating the validity of 
Eq. ||3J. The error bars for the concentration vs. time plot 
are of the size of the symbols and are not shown for clarity. 
The horizontal lines in the inset show the estimate of the equi- 
librium constant (thick solid line) along with the estimated 
statistical error (dashed lines) , while the circles are the values 
of the ratio Cg(t) /ca(£) computed from the curves in the main 
figure. 



is the interaction Hamiltonian that couples the reaction 
coordinate to the environment, with g controlling the 
coupling strength. The quartic coupling between the har- 
monic oscillators in fl e7W (P, Q) is intended to highlight 
the fact that Eq. J3Jl does not depend on the environ- 
ment having purely harmonic interactions, as in the case 
of exact treatments of similar particle-bath systems . 

Several parameters (a, b and <?), integration times t, 
temperatures, and numbers of oscillators have been used 
to confirm the validity of Eq. JSJl in the above model; in 
all cases, this equation has been verified to within the 
allowed statistical errors. Figure shows a particular 
set of results with a = 6, b = 1, g = 0.5, N = 10, 
and temperature T = 5, which were equilibrated with 
the Metropolis Monte Carlo method respecting the local 
equilibrium condition, and evolved in time with Verlet's 
Molecular Dynamics algorithm [Tij . For the experiments 
reported in Fig. |n = 100, 000 trajectories with time- 
step dt — 0.005 have been used for estimating the time- 
dependent concentrations, and the same number (n) of 
statistically independent samples was used for the esti- 
mate of the equilibrium constant. The choice of the pa- 
rameters a and b has been inspired by the relevant energy 
scales involved in the trans-gauche isomerization energy 
landscape of butane (see e.g. 0), and leads to a po- 
tential where species B (q > 0) is energetically favored. 
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Note that kT is comparable to the height of the poten- 
tial barrier separating the two wells, leading to a clear 
departure from the usual exponential behavior in first 
order chemical kinetics. The inset of Fig. ^ illustrates 
the fact that Eq. ® remains valid even in such cases. 

Finally, it should be noted that the independence of 
Eq. on the specific intermediate species of the reac- 
tion makes this result useful both as a predictive and 
as a diagnostic tool for complex reactions, such as pro- 
tein folding-unfolding. As a predictive tool, it allows one 
to directly infer the detailed kinetics of a reaction (i.e. 
its time-dependent concentration) in the direction oppo- 
site to that available from the experiment, provided of 
course K is known. This could be of particular inter- 
est in numerical simulations of protein folding-unfolding, 
which are computationally expensive and are often car- 
ried out in one direction only, typically in the unfolding 
direction so as to take advantage of the greatly reduced 
unfolding time-scales at high-temperatures and of the 
availability of well-characterized native states 0, H^ |. 
As already pointed out earlier in the paper, Eq. © im- 
plies that concentrations measured in opposite directions 
should have identical temporal profiles, and in particu- 
lar they should share the same non-exponential features 
(if any). This conclusion seems to be in disagreement 
with recent stopped- flow experiments with large (> 300 
residues) proteins m which unfolding was seen to 
be a simple exponential (two-state) kinetics, while re- 
folding presented appreciable non-exponential behavior 
consistent with the existence of at least two reaction in- 
termediates. This discrepancy can be explained by the 
asymmetric character of these experiments: In the un- 
folding direction, the initial population of proteins is in 
equilibrium with the buffer only, and the reaction evolves 
in a buffer-denaturant environment, while in the refold- 
ing direction the initial population is in equilibrium with 
a buffer-denaturant solvent, and evolves essentially under 
the effect of the buffer only. Thus, the lack of agreement 
with Eq. © in this context diagnoses that the introduc- 
tion (dilution) of denaturants to induce unfolding (re- 
folding) gives rise to reactions at fundamentally different 
conditions, which implies that these folding and unfold- 
ing experiments are not simply opposite directions of the 
same reaction. In this the experimental condi- 

tions of refolding are presumably closer to native condi- 
tions than those of unfolding (in the sense that in the for- 
mer the dynamics of the protein takes place practically in 
the absence of denaturants), one could in principle esti- 
mate the detailed unfolding kinetics in native conditions 
by feeding such refolding measurements and the equi- 
librium constant of the reaction into Eq. © ■ Neglecting 
the residual presence of denaturants during refolding, the 
quality of this estimate depends solely on how closely the 
population in the presence of denaturants approximates 
the unfolded subpopulation in native conditions. 

In summary, it was shown that the property of micro- 



scopic reversibility in chemical kinetics leads to a useful 
symmetry relation (Eq. J3J) between experiments hav- 
ing the same equilibrium constant but different initial 
conditions. Though it resembles the well-known detailed 
balance condition of chemical kinetics (Eq. J2J), these re- 
lations are different in scope and use: Whereas detailed 
balance is an inherently phenomcnological relation (inso- 
far as it involves rate constants) with applicability lim- 
ited to reactions well-characterized in terms of their el- 
ementary steps, the symmetry relation herein derived is 
independent of the phenomenology or of the intermedi- 
ate complexes that underlie the reaction, making it par- 
ticularly useful for complex reactions, such as protein 
folding-unfolding, as described above. Although not pre- 
sented here, analogous results can be easily derived for 
multi- molecular reactions. 
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